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Mathematical non-linear models for the boiling 
water reactor and the turbine sub-assemblies, studied 
earlier, have been modified and assembled to represent 
the Big Rock Point nuclear power plant in order to obtain 
better plant response, in the present work. Modifications 
intioduced are (i) the non-linear representation of the 
throttle valve (ii) the incorporation of the recirculating 
pianp dynamics in the reactor recircu].ation loop (iii) re- 
presentation of the voids by first order sweep model. To 
gain confidence in this modified plant model, tests like 
pump trip, control rod oscillation are simulated on the 
digital computer and are compared with the plant test 
results. Many other tests, whose plant test results were 
not available, are also simulated for better understanding 
of the plant response. It would be observed that the 
present model behaved favourably while comparing vjith 
the available test results. The thesis also reviews the 
work done in the field of modelling and optimal control 
techniques by the researchers and also presents recommen- 
dations for further work. 



CHAPTER 1 


INTRODUCTION 

1.1 IMPACT OF NUCLEAR POVfflR GENERATION 

The biggest problem of the word today is either 
to search new sources of energy or improve upon the already 
existing available sources of energy to cope up with the 
fast growing demand for energy. The sources of energy like 
coalj petroleum, geothermal, gas etc. may not be able to 
meet this growing demand mainly because of the fast dete- 
riorating world reserve for these fuels. This leads one 
to think of the nuclear energy as an alternative to these 
energy sources. Today, both the developed and developing 
countries are striving in a big way to master this source 
of energy. For instance it is expected that in U.S. in 
1985 the contribution by this source would be 37/^ of the 
total power against 0.5 percent in 1970* In India also 
two nuclear power stations are already in operation at 
Tarapur and Rajasthan and two more at Kalpakam and Narora 
would come up very soon. 

1.2 TYPES OF NUCLEAR POWER PLANTS 

The major constituents of a nuclear power 
reactor are the fuel and the moderator situated in the core 
where controlled fission takes place. The fission in the 
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fuel gives 2,5 fast neutrons at energies greater than 

235 " 

1 MeV for every neutron absorbed in the U nucleus 
(average). The energy of fast neutron is degraded by 
elastic scattering in the moderator from 1 MeV to 0,025 eV 
approximately. The neutrons \jhich are in thermal equili- 
brium with the assembly are called thermal neutrons (approx. 
0.025 eV) and are ready for further fissioning to maintain 
the chain cycle. This process is achieved in different ways 
which entirely depends on the type of fuel and moderator 
employed. It is known that natural uranium has smaller 
fission probabilities using light water as moderator. So 
enriched fuel has to be used to build critical assemblies 
while building nuclear reactor with light water as the 
moderator. These reactors are called as light water reac- 
tors which can be subgrouped as (i) boiling water reactor 
where coolant is at lov/ pressure and boiling is permitted 
and (ii) pressurised water reactors where coolant is at 
high pressure where bulk boiling is not permitted. Natural 
uranium can be used as fuel with heavy water as the modera- 
tor which absorbs very little amount of neutrons. This 
is the CANDU system. 

Gas cooled graphite moderq,ted reactors are 
popular in United Kingdom and advance gas cooled heavy 
water reactors concept is being developed at U.K, to 
achieve higher efficiency upto h-0% or so. These reactors 
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are basically thermal reactors since neutrons as thermal 
energies are used for fissioning and to maintain chain 
reaction. 

Another class of reactors is the fast reactor 
where the neutrons at higher energy (1 MeV) are used in 
maintaining the chain reaction. This precludes the use of 
moderators. The lia_uid coolant like sodium, potasium etc. 
are used to remove the energy from the reactors mainly due 
to their higher thermal conductivity property. In fast 
reactors source of the neutrons can be captured in fertile 
materials to produce the fissile material vrhich after pro- 
cessing can be used as fuel again. These are called as 
breeder reactors and have been developed by France and 
Russia. 

1.3 OBJECTIVE OF TBE PRESENT STUDY 

The objective of this thesis is to study and 
combine the models of the Big Rock Point boiliif^ water 
reactor and turbine sub-systems. Separately studied earlier. 
On this the complete plant model tests like R 0 Pump failure, 
control *od oscillation etc, are simulated to compare the 
response with the test results in an effort to validate 
the model. Since comparisons between the two are poor, 
modifications are introduced on this combined model to 
get a better match. Besides, other tests are also studied 
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to give an insight into the plant behaviour* 3he content 
of each Chapter is briefly discussed below. 

The Chapter 2 discusses the details of Big Rock 
Point boiling water reactor. The experimental results used 
in this thesis for comparison pertains to this reactor. 

The literature in the area of dynamics and control of 
nuclear reactors is surveyed, in great details in this 
chapter. 

Chapter 3 is devoted to the development of a 
mathematical model based on physical lai^rs and empirical 
relations. Previous work on this subject is briefly dis- 
cussed and modifications of this model are presented in 
detail . 

The dynamic simulation of the model given in 
Chapter 3 for various disturbances forms the’ basis of the 
Chapter k-. Comparison of the plant test responses with 
the model responses for some tests such as pump trip, 
control rod etc. is also given. Finally the chapter is 
concluded with the recommendations for further work. 



CHAPTER 2 


DESCEEPTION, CONDUCTED TESTS AND OTHER SIMULATION 
STUDIES FOR BRP POWER PLANT 

2.1 GENERAL FEATURES 
SITE AM) LOCATION 

The Big Rock Point nuclear power plant site is 
located in Charlevoix country, between the towns of 
Charlevoix and Petoskey, on the northern shore of 
Michigan’s lower peninsula. The site property includes 
600 acres of gently sloping wooded and cleared land on the 
shore of Lake Michigan at the western extremity of the 
southern shore of Little Traverse Bay. The site is 228 
miles NNW of Detroit and 262 miles NNE of Chicago. 

2.2 PLANT FEATURES 

The Big Rock Point nuclear power plant consists 
of a direct cycle, forced circifLation boiling water reac- 
tor, a power extraction system and associated service 
facilities. The principal structures include: 

i) a 130 ft. diameter spherical containment vessel, 
ii) a turbine building, 

til ) a structure housing water intake facilities, 
iv) a ^0 ft. stack waste storage vaults. 

The containment vessel -houses the reactor, recirculation 
piping, pumps, steam drum, fuel pool, and equipment for 
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removal of shut down heat. The turbine-generator and other 
conventional plant equipments are housed in a separate 
building. The containment is a spherical steel vessel of 
130 ft. dia, its prime purpose being to prevent any haimfuL 
spread of radioactive material to the environment in the 
event of an accident in the reactor system. It also serves 
as a weather-proof housing, 

2.3 PLANT DESCRIPTION 

For detailed description and design criteria one 
shoiJld refer to references |l “ ^ maj'or components (Fig. 2*1.) 

are. briefly described below. ■ • 

REACTOR VESSEL 

The vessel is cylindrically shaped with full 
hemispherical top and bottcan as shown in the Figure 2.2. 

It is 30 ft. long, 106 inch in diameter with wall thick- 
ness 5*25 Inch and cladding thickness 5/32 inch. It is 
designed for 1700 psia pressure, 650 °F temperature, where- 
as normal operating conditions co-uld be 1500 psia, 593 

REACTOR VESSEL PBNETPATIONS s 

Penetrations Number Diameter Location 

(inches') 

Coolant water inlets 2 20 Bottom 

Steam-water mixture outlets 6 1^ Shell 

Outlet to shutdown heat exchanger 1 6 aiell 
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Access ports in top head 3 
Control-rod drive penetrations 32 
Liquid poison inlet 1 
Emergency core spray inlet 1 
Vessel vents 2 
In-core flux monitors penetrations 8 
Instrument nozzles 4 
Seal leak monitor 1 


10 Head 
4 Bottom 
3 Bottom 
3 Shell 
3 Head 

2 Bottom 

3 ^ell 

1/2 Head 

flange 


FUEL AND COEE DESCEIPTION 

Each fuel bundle contains 144 fuel rods in square 
array, 132 being standard rods and 12 being special rods 
used to minimize the local hot spot factors. Ihe fuel 
channels surrounding each fuel bundle provide guide surfaces 
for the control rods and allow orificing of the flow between 
channels. The channels are supported by support tubes. 

There are 32 control rods cruiciform in nature (11 inch 
wide X 5/l8 inch thick) consisting of 116 poison tubes 
containing B^ G powder. When control rod movement fails, 
scram is initiated by pilot valves to make the reactor 
sub-critical till it cools down, the liquid poison (850 
gallons) being located above the steam drum and consist- 
ing of sodium pentaborate solution i'S &2 0^^ 9 H 2 O) 

obtained by dissolving boric acid and bo3?ax at 212 ^F, 
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BIOLOGICAL SHIELDING 

To protect plant personnel from radiation 
hazards due to neutron and Z - radiation, a hiological. 
shield is provided based on the international dose accep- 
tance level of 5 rem per year. It is made up cf heavy 
concrete and is provided with necessary cooling. 

STEAM DRUM 

Hie steam drum is situated high up in the enclo- 
sure. Its main function is to separate steam from the 
steam-water mixture from the reactor and provide exit 
steam of 99 *9/^ quality. With its storage capacity of 500 
cubic feet, surges of water level and pressure between the 
reactor vessel and the drum are easily accommodated. It 
also assures NPSH (net positive suction head) for the 
recirculating pumps. Mater level in the drum is 65 feet 
above the pump suction level and provides adequate natural 
circulation driving heads, when the piMps are free to 
rotate but not being driven by the electric motor coupled 
to it. It will be seen ( Sec 1 .1) that the reactor under 
such condition returns to ^0^ of normal power*. The drum 
has an axial length of 4-0 ft, with ID of 78 inches. 
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RECIECI3LATIMG PUIPS 

There are two vertically mounted recirculating 
pumps operating in parallel, each capable of returning 
17000 gpm of water to the core. The pumps are single- 
stage, centrifugal, double-suction with an overhung impeller ♦ 
The pump driver is a vertical, drip-proof induction motor 
of conventional design, rated for hOO HP at I 6 OO K’M, 2300 
volts . 

PRIMARY LOOP PIPING AND VALVES 

The primary loop piping interconnects all the 
major components of the nuclear steam supply system. Six 
inch risers carry the steam-water mixture from the 
reactor vessel to the steam drum. After the steam is sepa- 
rated, the recirculating water flows out of the drum down 
four 17 in. downcomers. The downcomers, in two groups of 
two each, join together into two ^ in, pump suction 
headers. The two pump discharges are each 20 in. and return 
to the bottom of the reactor vessel. A system of vibration 
dampers are fastened to the piping to reduce possible 
vibration stresses. An open bypass between the suctions 
of the two pumps will keep water flowing in the four down- 
comers even if one pump fails essentially, overriding the 
possibility of tilting of the drum due to unequal expan- 
sion in the two sets of downcomer piping. Each pump has 
24 in, gate valve at suction and 20 in. gate valve at 
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discharge end to keep out of service if leakage becomes 
excessive. In addition a 20 in. butterfly valve at the 
discharge end regulates flow during research and develop- 
ment programmes. When the valves are fully closed, scram 
is initiated. Six safety valves are located on the main 
steam drum meant for operation beyond I870 psia pressure 
in the drum. 

TUEBIliE GENERATOR EQUIPMENT 

The turbine is a 3600 rpm, tandem compound, 

double flow, condensing unit, directly connected to a 

hydrogen cooled generator v/hich in turn is connected 

through a reduction gear to an air-cooled exciter. Three 

points of extraction for feed water heating are provided. 

The turbine is rated at 5^,500 KW at 1000 psig. 0 degree 

-] 

final supeiheat and 3 2 mercury absolute exhaust 

pressure, with 3 % make-up a3_lo\fance and three feedwater 
heaters in service. The turbine is capable of operating 
continuously at lif 50 psig, 0 degree final supeiheat and 
2 inches Hg back pressure with a maximum expected output 
of 75>000 KW. The 13j800 volts, wye-connected generator 
is rated at 7 O 5588 KVA, 0,85 power factor, 0,8 short ' 
circuit ratio at 30 psig hydrogen cooling pressure. 

Besides conventional design criteria, all modi- 
fications necessary due to use of saturated steam from 



the BWR are incorporated in the turbine design. Particular 
attention is paid to eliminate pockets or crevices in which 
radioactive material may lodge. Each turbine stage is 
drained either internally or externally. The turbine is 
provided with moisture removal buckets ahead of each extrac- 
tion point. In addition, two external moisture seperators 
are provided in the cross-over between the high pressure 
and low pressure sections. 

MAIN TURBINE CONDENSER 

The main turbine condenser performs the following 

functions s 

a) condenses the steam exhausted from the turbine to 
obtain the desired heat utilization and vacuum, 

b) de-aerates the condensate and water from heater 
drains and other returns, 

c) serves as a heat-sinli for excess reactor steam 
dumped through the turbine by-pass valve, and 

d) detains condensate in the hot well to permit 
decay of short lived radioactivity. 

The condenser is made-up of steel and is a horizontal, 
single pass, divided water-box de-aerating type unit of 
conventional construction. It provides 27 j 500 square 
feet of effective condensing surface area. A 6000 gallon 
oversized, baffled, storage- type hot well is provided to 
allow decay of short-lived radioactivity. 
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CONDENSATE AND FEEI> WATER SYSTS4 AUXILIARIES 

a. EXTRACTION DRAINS AlID VENTS 

Extraction steam for feedwater heating is taken 
from three points off the turbine to the respective heaters. 
The two higher pressure extraction lines are provided with 
automatic bleeder trip valves to protect the turbine from 
flooding in the event of a heater tube break or overspeed 
from steam flashing out of the heater after a turbine trip. 

Water collected from the turbine moisture removal 
stages is piped to the drain cooling section of either the 
high, intermediate, or low pressure heaters. Heater drains 
are cascaded to the condenser where they are de-aerated 
and merge to the condensate flow. 

b. CONDENSATE PUMPS 

Two half -capacity, vertical, multistage centri- 
fugal pump* pumps the condensate from the hot well through 
the condensate systeiii to the suction of the reactor feed 
pumps. 

c. FEEDWATER PUMPS 

Two feedwater pumps take suction directiy.i from 
the condensate system, ciLschai^e feedwater through the 
high pressure heater and through a common header to the - 
reactor steam drum* These are horizontal, multistage, 
centrifugal pumps. 
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d. FEEDWATER HEATERS 

Three feedwater heaters are located in the conden- 
sate circuit. The flow pressure and intennediate pressure 
heaters are of the horizontal mounted, U - tube type with 
removal tube bundles, integral drain coolers, and bolted 
head covers. 

TURBINE AND BIPASS CONTROL SYSTEM 

The turbine obe 37 -s two modes of control? 

a. Initial pressure regulation (IPR, Base-load operation) 

b. Speed control (House-unit operation). 

Normally, the Turbo -Generator set is base loaded and to 
avoid swings with accompanying transients being felt by the 
reactor, the steam main line pressure is maintained cons- 
tant by IPR, This essentially, controls valve position 
without regard to the generator or transmission line system 
loads. But in the event of excessive speed build-up, IPR 
is overriden. This excessive build-up noimally happens . 
when the I38 W transmission line trips (due to some trans- 
mission faults) and load is rejected on the generator 
actuating turbine speed-up. 

AUXILIARY COOLING WATER SYSTEMS 

The service water system is an open system in 
which strained water is supplied from pumps in the intake 
structure and returned to the lake along with the discharge 



from the circulating water system^ essentially, to cool the 
following equipments 

Generator hydrogen coolers, 

Turhine tube oil coolers, 

Feed pump hearings and oil coolers, 

Air compressor, after-coolers and jackets, 
Miscellaneous sample coolers, 

Airconditioning system, 

Reactor enclosure air coolers, 

' ^ace heating and cooling equipment, 

Reactor Cooling water heat exchangers. 

REACTOR COOLING vaiER SISTM 

This is a closed intermediate cooling loop utili- 
zing demineralised water to cool the following equipments s 

Reactor shield cooling panels, 

Reactor clean-up Non- re gene native heat exchanger, 
Miscellaneous sample coolers, 

Reactor recirculating pump coolers. 

2.^+ DESCRIPTION OF TESTS CONDUCTED ON TEE PLANT 

These tests are extensively documented in 
references [4, 5] and are summarized below. 
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CONTROL ROD OSCILLATION TESTS 

The aim was to measure the closed loop, reactivity 
-to-neutron-flux transfer function of the primary hydraulic 
loop for various conditions of flow, power and pressure* 

One of the standard control rods was made to oscilla.te in 
the frequency range 0.2 - 5*0 cps, while keeping reactor 
steam pressure controller, steam ~ drum level controller 
and the movement of control rods (other than the oscillating 
rod) on manual operation to avoid unwanted variation of 
parameters due to controller cycling. 

NOISE TESTS 

The aim was to collect data of the noise in 
steady state flux and flow signals after performing rod 
oscillation test. These tests were of long duration (^5 
minutes) following immediately after CRD test. Feedwater 
level controller and reactor steam pressure controller 
were put on manual operation to eliminate noise due to 
controller cycling effects. 

PRESSURE OSCILLATION TESTS 

These were performed to detemine the pressure- 
to-neutron-power transfer function of the reactor primary 
hydraulic loop. The main steam valve bypassing the tur- 
bine was oscillated over the frequency range 0,02 - 
3.0 cps. 
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PEESSCJRE TRANSIENT TESTS 

These tests were performed tefore control rod 
oscillation tests. Ihe drum level controller was put on 
manual. Reactor pressure was controlled ty first opening 
the turbine bypass valve by 3 % and then changing this set 
point rapidly. The pressure, neutron flux, steam flow and 
feedwater flow were continuously monitored during tlie test, 

PIMP TRIP TEST 

The tests discussed so far involved small order 
perturbation, while the pump trip test and the otiiers that 
are discussed below involved large-order changes. 

For various operating values of pressure, tempe- 
rature and neutron po\'^er,the po’-'er to both recirculation 
pumps was tripped and signals were continuously recorded 
for reactor pressure, neutron power, average heat flmc, 
core inlet velocity and core inlet temperature. There 
were several variations of the pump trip test, e,g., 
tripping of a single pump, restart of the pumps after the 
trip etc. 

FLOW CONTROL TEST 

With the reactor operating at rated flow and 
at a certain power level, the recirculation flow was 
decreased in steps to 90, 80, 70, 60 and 50 percent of 
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the rated flow without the control rod movements. The data 
so collected helped in thorough understanding of core 
performance, 

LOAD FOLLO’rflNG TEST 

Mith the plant operating at specified operating 
point, the turbine was put on governor control and feed- 
water on automatic operation. Electrical load was rejected 
from the turbine at rates of 1, 2, 3, 5 MWt per minute. 

For each rate, pressure regulation was obtained by manual 
throttling of the butterfly valves in the recircula-tion 
loop. During this test, control rods ^^ere held at fixed 
positions. 

NATURAL CIRCULATION TEST 

This test was conducted for the 84 bundle core 
configuration only at 1350 psia. The plant was Initially 
operating at 75 MUe at rated two-loop forced circulation. 
Die control rods were inserted in a certain specified 
pattern. After the reactor steadied out, the pumps were 
tripped till power and flow coasted down to their natural 
circulation levels. Thereupon, power level increases 
were affected in accordance with a specified programme 
of control rod movements. For each new control-rod 
pattern, data was obtained to calculate the reactor power 
from the heat balance equations. Prior to each power 
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increases, a comparison was made "between the actual 
natural circulation flow and power conditions and the 
licensed heat transfer limit, Vftien the comparison indica- 
ted that it was permissible, further power increases were 
made untill the heat transfer limit ^^as nearly reached. 

2.5 SAFETY STUDIES 

DIUA14IC STABILITY OF THE EEACTOE 

Some of the low pressure ( — ■ 500 psi) experimen- 
tal BWR' s exhibit power oscillations due essentially to 
mutual feedback between hydrodynamic flow oscillations and 
the effect of voids on reactor power. The design criteria 
for BRP which ensure stable operation include high system 
pressure, forced circulation, a conservative void reacti- 
vity coefficient, a long fuel heat- transfer time constant 
and the use of a pressure regulator to control system 
pressure. The analytical studies reported in GEAP-3795 L6J 
based directly on the physical concepts of momentiaa inter- 
change, energy conservation and mass continuity, also do 
rule out the possibility of any such instabilities. 

Usual unsafe operating conditions for the plant 
would be short reactor period, high neutron flux, hi^ 
reactor pressure, low water level in the primary steam 
drum, simultaneous closure of turbine stop valves and 
primary by-pass valves, closure of backup primary steam 
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sphere isolations valves, high enclosure pressure, low 
condenser vacuum, low water level in the reactor vessel, 
sim-ultaneous closure of the recirculation water line valves, 
and high level in the scram dump tank. At the onset of such 
potentially unsafe conditions, the reactor is automatically 
scrammed and cooling initiated through the emergenci^ con- 
denser. 

The experimental tests discussed in Sec. 2,4 
serve to varify the inherent stability of the plant. 
However, one should consider the safety requirements in 
the event of equipment malfunctioning, operator error or 
some combination of these. 

SAFETY DESIGN AGAINST EQUIPMENT MALFUNCTION 

The design requirements for safety necessitate 
the following i 

i) use of equipment with the highest practical 
reliability against failure, 

ii) fail-safe design, i.e. any failure will result 
in action toward safe direction, 

iii) safety througli redundant devices (mostly 

triplicate system en^^loying 2/3 redundancy, 
is in voguejk 

The possible important equipment failures coifLdbei 
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CONTROL ROD MALFUNCTION 

The following may help in its initiation: 

(i) failure of normal drive power, (ii) failure of emer- 
gency drive power, (iii) seperation of rod from drive. 

In all such cases, reactor is shut down through scrams. 

LOSS OF ELECTRICAL LOAD 

A sudden loss of electrical load will cause a 
partial closure of the turbine control valve. The bypass 
steam valve will open automatically to hold reactor system 
pressure constant by providing a path to the condenser for 
the turbine- rejected steam. Thus, there will be no unde- 
sirable flux transient from the load change. The control 
rods would be manually repositioned to hold reduced power. 
In the event of bypass valve failure, or its slow operation, 
the reactor will be scrammed through the high neutron flux 
monitor or through a high pressure signal bringing up 
emergency condenser action, in case the flux monitor fails. 
In the event of canplete load rejection the reactor is 
scrammed if both the turbine and bypass valve remains 
closed. 

LOSS OF CONDENSER VACUUM 

If the condenser vacuum is lost, the main reac- 
tor heat sink is lost. If no action were taken by an 
automatic protection system the condenser pre sure would 
increase to a point where the condenser rupture diaphragm 
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would fail to ease off this built pressure. This possibi- 
lity is eliminated by sending a scram signal to the reactor 
through a vacuum sensing device. Also the turbine is 
tripped. Othen^rise, if the device is inoperative, closure 
of the bypass and stop valves is demanded and this intro- 
duces scram directly or indirectly through high pressure. 

If the cause were loss of coolant flow due to power failure, 
the reactor receives a shutdown signal through an overall 
energised safety system. 

FUEL CLADDING FAILURE 

■ An off-gas hold-up system provides protection 
when fuel cladding fails. Hie system is designed for 
manual closure of the gas holder isolation valve at a noble 
gas concentration that could be reasonably expected to 
deliver the maximum permissible off site dose if the opera- 
tion were to continue inabated for a year. Automatic shut- 
down will occur at a greater gas release rate in time to 
keep personnel exposures within acceptable limits. 

LOSS OF FEEDWATER HEATERS 

Sudden loss of all the feedwater heaters would 
cause an immediate but smooth rise in flux. lAfithin a few 
minutes the flux would be high enough to initiate scram, 
and no overshoot or damage to fuel would result. 
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LOSS OF FEEDWATER 

Loss of feedwater will I'esixLt in gradual lowering 
of tile water level in the steam dTum and, if continuedj 
will automatically initiate reactor shut down. Main 
condenser cools the reactor adequately, 

SaiFETT DESIGN AGAINST . OPERATOR ERROR 

The probability of accidents being initiated 
by operator mistakes is very less, considering the level 
of training administered to reactor operators. Nevertheless, 
safety is ensured by providing detailed operating and main- 
tenance procedures, and by special design measures for 
control of certain operations* For example, the reactor is 
never started up unless all safety systems have been properly 
actuated. Further, only one control rod can be withdrawn 
at a time and that too at a limited rate. In PSAR C2J a 
number of hypothetical situations are discussed which 
pessimistically combine operator error with equipment fai- 
lures to result in 

a) start up accident, 

b) cold water accident, 

c) fuel loading accident, 

d) closure of steam line back-up isolation valve, 

■ e) failure to replenish cooling water in emergency 


condenser. 
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It is shown how such accidents are rendered impossible by 
the various safety systems in BEP. 

SIMULATED ACCIDENTS 

WASH - 1270 and WASH - 1400 describe a 
variety of simulated maximum credible accidents for nuclear 
power plants, resulting from 

a) primary system rupture, 

b) sphere pressure reduction, 

c) leakage from enclosure, 

d) large LOCA and small LOCA. 

It is shoi-m in these reports that such occurren- 
ces have extremely low probabilities. References ["93 lol 
on the maintenance and operational record of BEP are proof 
of the excellent safe performance of the power plant to 
date. 

2.6 REVIEW OF EARLIER LWR SIMULATION km CONTROL STUDIES 

A non-linear model of the reactor and turbine for 
load following studies of BRP has been presented in IBM 
research report fl1 3 GIRIJASHANKAR p2, 1,3 3 1^”] consi- 
dered ttie simplified 8 - order linear version of the 
nuclear reactor to simiilate the small order disturbances 
separately and then by regular approach studied feedback 
control laws that drove the systan to steady state. 



GIEIJASHAMAR also studied an 11 - order linear model 

of the nuclear turbine for small order disturbances. 

LINROED [l 6 j| made a GEG BWR model for transient studies. 
GLAASSEN and ECKERT fiyjj in the suggested direction of 
WASH - 1270 simulated tests without scram to see the 
plant response of a BWR more rigorously. Earlier BECKER 
{j:8, 19jstudied bumout conditions in channels both experi- 
mentally and analytically. BECKER et al. [aol also pointed 
out conditions for hydrodynamic instability, LEONARD, SUN, 
MUNTHE ANDERSON and DIX [" 21 J developed a GEC model for 
analytical LOCA. studies- for reactor safety evaluation. 

MEHTA [^223 studied non-linear simulation of an organic- 
cooled Canadian CANDU reactor. CHOU (”233 considered a 
detailed model of Pickering G,S. 'A* for load following 
studies., KERLIN [ 2^3 simulated H.B., ROBINSON BWR nuclear 
reactor for reactivity and steam valve perturbations. 

ATRAY and SHAH £253 modeled a complete prototype power 
plant for analytical studies. 

Such studies of modelling and simulation nece- 
ssarily validate the model before conducting optimal 
control studies.' FRDGNER £263 has highlighted the impor- 
tance of optimal control studies and in his Ph.D. thesis 
[^1 has given a useful reviex-j. He himself applied feed- 
forward technique on a 13 order linear model derived 

from 25 order non-linear differential equations using 



Kalman filter technique for non estimable state variables 
to implement Direct Digital ContiV)! on BWR. The algorithm 
was tested on estimation of void reactivity coefficient 
under different conditions by RA.0 [28J. FROGNER [29^ has 
surveyed optimal control methods and discussed [30I iden- 
tification techniques in the presence of plant noise. 
GRUMBAGH [_3lJ stressed the simultaneous need of software 
and hardware control structures less sensitive to component 
failure in control and instrumental failure, rie built up 
the conceptual design for such techniques. Modern concepts 
are more or less based on Kalman Filter techniques and most 
of the work simplified to date is nucleated around these. 
Thus, GODBOLE studied the application of Kalman 

filter technique in nuclear reactors to obtain the non- 
measurable variables. He later [ssHi applied Kalman filter 
technique in the presence of white process contaminated 
with noise, 

VENERDUS and BULLOCK used Kalman filter 

technique to estimate time dependent reactivity. But this 
technique increased the computational burden, and SAMANT 
and SOREUSON devised a low-order filter to save 

computer time. Little attended observers were studied 
LUENBERGER [36, 37^ . GODBOLE |_383 designed a new pro- 
cedure which (a) took care of linear feedback laws without 
sacrdficing observer poles, (b) had a true deadbeat 
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observer which tracks the acbual states iii observability 
under time steps. 

Presently, research interests are to apply 
process computers for the optimal on-line control of 
nuclear power plants. LUI'TDE [39, 4 -oJ has discussed the 
development and use of such hardware techniques. 

BLOMSNES ^2^ considered the power behind linear state 

variable feedback control wliile considering spatial core 
distributions. . GRUliBACH [j+sJ. proposed the computer usage 
at Hal den for controlling the spatial distribution of 
neutron flux. GRIMBACH [VfJ analysed plant disturbances 
through process computers. GRUMBACH and BLOMSNES jj+53 
presented practical applications at Halden, BLOMSNES and 
NETLAND [ifbj proposed the use of a simulator and trained 
operators for computer based reactor control. JOSEESSON 
[ 4-7 J in a cooperative effort between Sweden and Norway, 
showed the successful use of advanced control methods at 
Studsvik through the simulator STUDS-1 where load-follow 
control studies were conducted on the Nordic electric 
power grid. BJORLE et al. 4^, BLOMSNES et al. '^6^ 

GRUMBACH |_51 "j have reported the recent application of 
optimal control laws to HB¥R at the OECD Halden Proiect, 
highlighting the advantage of DDC control over conven- 


tional methods. 
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OBJECTIVE OF THE PEESEHT WORIv 

Most of the work reported till date ignores the 
turbine dynamics or represents it in a simplified foim, 
as M/llJKIN and SHINOHARa £72^ This is good enough when 
boiler transients are initiated by weak disturbances and 
persist only for small duration. But for large order 
disturbances lasting for minutes, the turbine dynamics 
cannot be ignored or represented by its associated time 
constants by simple equations. The present effort is 
concentrated to develop the analytical non-linear model, 
introducing suitable modifications, for the BWR power 
plant suitable for large order disturbances. 

This kind of models could be useful in ' 

a) building simulator employing real-time process computer 
to train the operator to have good feel for both normal 
and abnormal opera.tion and the plant behaviour, 

b) to conduct optimal control studies, 

c) to check the plant behaviour under anticipated tran- 
sients without scram. 

The present study discusses the nonlinear model 
of the turbine and the reactor subsystems separately. 

Then on the combined model which represents tiie boiler 
and turbine, the response for the lai^e order transients 
are studied. The BBP boiling water reactor plant data 
will be used to validate the model. 
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CONCLUSIOIS 

The important components of BHP boiling \^ater 
reactor and the tests conducted on it were discussed* 

Some, of these will be simulated and discussed in Chapter h » 
A useful literature survey giving information on modelling 
and control technique was listed. Finally the object of 
this work was outlined. 
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CHAPTER 3 


NON - LINEAR MODEL STUDIES FOR BRP POWER PLAITT 

3.1 GENEPAL CHARACTERISTICS ' 

We dealt v/ith the plant description in Chapter 2 

for the BRP boiling water reactor. The salient features 

of its operation are reviewed. For details one should 

refer to the texts such as Refs. [53 - 

235 

core, the fission process of the U nucleus occurs in the 

fuel releasing 200 MeV per fission. The major part of the 

heat generation is via the slomng down of the fission 

fragments within the fuel, A small fraction of the heat is 

contributed by neutron and gamma heating of the coolant/ 

moderator, structural materials etc. The neutrinos escape 
' » - 

the systenf with a fraction of irrecoverable amount of the 
total energy i.e. 200 MeV. The transport of heat from fuel 
to canning is by conduction and from canning to the coolant 
(light water) is by convection. This heating of the coolant 
converts the single phase coolant into a two phase mixture 
in the core. Steam in this two phase mixture is separated 
by steam separators in the drum when the two phase mixture 
advances into the drum. From here the steam is flowed to 
flow into the turbine as shown in Fig. 3.1, Once the 
steam is seperated in the drum the remaining saturated water 
in the recirculation flow mixes with the feedwater 


55 ]. In the BWR 
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(subcooled) in the drum and then flows down the downcomer 
as subcooled water to enter the inlet of the lower plenm 
chamber and then the core throiogh the recirculating pumps* 
20 % of the recirculating water at inlet cools the control 
rods while the rest 80^ cools the fuel elements. 

Steam from the drum enters the steam chest of the 
turbine through the throttle valve and expands in the HP 
unitj a part of this steam is bled out at exit. At this 
point the moisture separators tra.p the moisture to increase 
the qualit!” and returning the collected moisture to heater 
No, 2 while the main steam enters reheater. The reheater 
superheats the steam which then expands through the IP <& LP 
■units. A part of this steam is bled out and the rest is 
passed to the condenser. The condenser also receives the 
main steam through by-pass valve in emergency conditions 
and in conditions where tests like valve rod oscillations 
are conducted. The condenser operates at uniform pressure 
(ass-umption in our study) . The water, coming out of conden- 
ser is heated in reheaters No. 1 and 2 by bled steam out 
of IP & LP, and HP unit. This heated feedwater returns to 
the drum where it takes 6 seconds to mix-up with the 
saturated water of the recirculation flow. 
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3.2 BASIC MODEL EELATIOWSHIP 

These physical processes are translated to mathe- 
matical language (Sec. 3.3) through four fundamental laws 
viz .s 

a. CONSEHTATION OF ENERGY 

Energy is neither crea.ted nor destroyed, if spent 
in one form it reappears in other forms. This essentially 
implies that the time rate of change of energy in a control 
volume is given by the balance between its generation rate 
within the control volume, the entrance rate from the 
surroundings and the exit rate to the surroundings. 

b. CONSERVATION OF MOMENTA 

The net force of an inertial layer is equal to 
its rate of change of momentum. 

c . CONSERVATION OF MASS 

For each component the continuity equation is 
preserved, i.e. the rate of mass increase in a control 
volume is equal to mass flow rate at entrance less that at 
exit. 

d. EQUATION OF STATE 

For any thermodynamic cycle which consists of 
thermodynamic process of only PdV t3q)e of work mode 
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there are only two imdependent variables such as pressure 
and enthalpy. Any other variable can be represented in 
term of these two variables, Reynolds f58'"| , This is the 
State Postulate Law. 

3.3 PREVIOUS MODELS 

3-. 3.1 THE REACTOR 

GIRIIASHANKAR jj2, 13_j essentially considered the 
IBM non-linear model jjO? linearised it for the control 
studies of the reactor. The basic development x^ras in accor- 
dance with Sec. 3*2, with the following assumptions : 

(i) The model does not calculate minimum critical heat 
flux ratio (MGHFR) or the maximum fuel temperature 
in coolant channel, local velocity changes and local 
’toid changes etc. Fast transients were not included. 

(ii) The entire neutron population responds as one energy 
group. The spatial variations of the neutron flux 
are ignored, and a single delayed-neutron group 
being considered in representing the core kinetics. 
Only the reactivity feedback due to changes in fuel 
temperature and void changes are considered and no 
control rod sequential pattern insertion is studied. 

(iii) The core average void fraction is a constant frac- 
tion of the core-exit void fraction. 
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(iv) The tempeiature distribution of the fuel is indepen- 
dent of axial and angular space variables. Fuel 
density and specific heat are constant across the 
fuel element. Storage of heat in fuel-clad gap or 
heat generation in the clad is negligible. 

(v) In the hydro djmamic study, the variables like core 
steam quality and void fraction vary instantaneously 
with the changes in coolant flow, coolant inlet tem- 
perature and heat flow to the coolant. Pressure in 
the steam drum is spatially invariant. The steam and 
water phases are always in thermodynamic equilibri-um. 

(vi) In the derivation of the momentum equation, compre- 
ssibility effects in the loop are ignored. Slip, 
the ratio of average fluid velocity to the average 
steam velocity at core exit is assumed independent 
of void fraction and mass flow rate, 

(vii) The recirculation flow liquid is subcooled, incompre- 
ssible and one- dimensional flow in the loop is 
assumed, i.e. the velocity heads ir the bulk water, 
downcomer and lower plenum are in the direction of 
flow. The core leakage flow remains at the core 
inlet enthalpy until mixing occurs at the core 
outlet . 
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(vlii) The flashing process occurs instantaneousl 3 '' after 
pressure changes in the system when the two phases 
are in thermal equilibrium. 

(ix) A transport delay of 9 seconds between the drum 

and the inlet plenum chamber is included disregard- 
ing delay due to mixing in the steam drum. 

V/ith these assumption, the model so developed is 
represented as follows; 
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fg = 2.25 X 10“^ P ibm/ft^ (3.^2) 

h^g = 6.85 X 10^ - 2.08 (P - 1000) x 10^ W Sec/lbm 

(3.^3) 

hj = 5.72 X 10^ + 1.81 (P -■ 1000) X 10 W Sec/lbm 

(3.Mf) 

Equations 3*^0 to 3.^i-^ 3-ne empirical relations 
around 1000 psia deduced from Mollier charts and steam 
tables. 


Eq^uations 3-.16 aii.d3*17 are due to Becker’s expre- 
rimental correlation 58j where he observed that ’A’ is 
independent of mass flow rate. This, however is not strictly 
true, but for the present study we will use Ms foimtila for 
the two-phase friction multiplier. 


3.3.2 TUEBINE MODEL 

GIRIJASHANKAR jl53 compared the validity of non- 
linear and linear turbine models which were developed 
based on IBM report [l{l* 

The development assumed the followings 


(i) 

steam inertia and 

kinetic energy changes are 


ignored, 


(ii) 

a perfect-mixing, 

adiabatic process prevails 


the steam chest, 
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(lii) th./e throttling is a pure isenthalpic process 

allowing no change of enthalpy across the valve. 

For a given valve stem position, the flow tiirou^ 
the valve is proportional to the upstream pressure 
of the valve, 

(iv) bleed flow to the heaters for regeneration purposes 
is at reheater pressure and so is a fraction of 
instantenous flow through the turbine. iJeither 
representation of theimodynamic path nor calcula- 
tions for individual stages are attempted, 

(v) the overall efficiency of the HP, LP and IP turbine 
is modified to yield a luliped efficiency. Correc- 
tion factors are used to take into account 
efficiency changes due to stage leaks, root and 
trip losses, rotational losses and reheat effects, 

(vi) the dynamics associated with the inteimediate and 
low pressure turbine are secondary to the dominant 
dynamics of the reheater. Thus IP & LP are repre- 
sented in a lumped single unit, 

(vii) the condenser pressure is assumed to be invariant, 

(viii) full flow is used to produce net V7ork through the 
HP turbine while one-half of tiie bleed flow- 
produces torque in the IP and LP turbines. This 
is because bleeding is done at the end of HP 
turbine while in the IP & LP turbines, it is 
carried out at various points. 



With these assmptions the follo\d.ng equations resiilt 


K3ZZLE BOWL DYNMIGS 


(W^ hg - hg) J pg 

pc\ EF 
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dt I 1 - k1 J 
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(3.^5) 
( 3 .^6 ) 


S ^1 ^ 


(3.^7) 
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(3.^9) 


- Fc ^G “ 


(3.70) 
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(3.51) 
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(3.53) 


MOISTURE TRAP 
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(3.58) 

( 3 . 59 ) 


Cb) REHEATER 'STEAl'i 
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( 3 . 62 ) 

(3.63) 
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%2P 

%2P 

OUTPUT EQUATIONS 



%P 

= 

Kp ^ijjp 

*** ^2^ 

(3.71) 

^LP 

r: 

% "*1lp ^3 

(hjj - h^) 

(3.72) 

"^HP 

=: 

^.CHiP HhP 


(3.73) 

>lLP 

= 

"^IcFLP ^LP 


(3.7^) 

\lp 

= 

^ HP “ hg 

-dip 

f _ 0.86 
^2 

(3.75) 

^GPHP 


\ 

(3.76) 
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"'ICFLP 
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( 3 . 77 ) 


(w' / Wn)(1 

- Wt ,) 



= 0.02 

TORQUE = Tgp + Tpp 
POWER = TORJ^UE x -Cl-. 

THE PARAJ-iETERS IN EMPIRICAL FOB! 
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(3.78) 

(3.79) 

( 3 . 80 ) 



0.018^ 


(3.81) 

“ 

2.288 - 0.0166 (P^ - 200) 


( 3 . 82 ) 

= I 

JI90 - 0.0256 (P - 1000 ) I 1 

.055 X 10^ 

( 3 . 83 ) 

hp = ! 

[355.6 + O.i+if (Pp - 200 )”j 1 

.055 X 10 ^ 

( 3 . 84 ) 

hfg = I 

[343 - 0.4 (Pp - 200 )[]i .055 

X 10 ^ 

( 3 . 85 ) 


hp + hpg 


( 3 . 86 ) 


Hp = (1067 + .37 (Pj^ - 200) - 0.0011 (Pp - 200)^ 

- 0.1 (Pq - 1000) 1.055 X 10^ (3.87) 

= 69.0 X 1.055 X 10^ (3.88) 



( 3 . 89 ) 

(3.90) 


3, If MODIFICATIONS INTRODUCED 

The non-linear models for sut-assemblies viz. 
reactor and turbine, presented above in Sec. 3.3.1? 

Sec. 3 . 3.2 were linearised around an operating point and 
validated for small order transients in earlier studies 
Refs. 11, 12, 13 , 15 . ®ie present interest is to 



validate and study the non-linear models for the entire 
assembly for large order transients, viz., RG pimp failure 
and turbine trip etc. In order to do this, the two sub- 
assemblies are coupled together and the simiiLation is 
carried out for the RC pump trip failure tgst. It is found 
that the response for neutron power is very sudden and the 
absolute values of various variables are far from satis- 
factory (Fig. 3.2). So it is felt that the model needed 
modifications and -tiie following describes the effort in 
this direction. These modifications are (a) better throttle 
valve representation (b) incorporation of the pump dynamics 
and (c) inclusion of sweep effects of voids. This study 
was presented in a recent paper by GIRIJASHANKAR, CHAVJLA 
and MALHOTEA where the improved effects were exclu- 

sively shown. These improvements are described as follows. 

3.1-K1 THROTTLE VALVE 

Isenthalpic throttle valve connecting ttie reactor 
to HP turbine was represented by a linear relationship 
between steam flow, area of flow and upstream pressure in 
Sec. 3.3.2 by Eqn. 3.^7. This is valid only if the down- 

m 

stream pressure is kept constant and the valve movemsnt 
is small in the neighbourhood of an operating point. 
However, in the large transient study, since the reactor 
feeds steam to the turbine the downstream pressure is 
likely to change and the valve movement is large enough 
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to introduce nonlinearities in the flow, pressure and valTe 
movement characteristics. Thus, the steam flow through 
the valve is n nonlinear function of the valve— stem posi- 
tion and difference of upstream and downstream pressure 
rather than upstream pressure. This new relationsliip is 
used while repre'senting the valve, in accordance to IBM 
report, 1972, as shown below 

N 

* C 7^ - y ) » (pg - Pp) (3.91) 

The constant is determined under the initialisation 

conditions, KEARTON |^55 | • I'he value of N is selected 
such that at rated flow, a ten percent closing of the valve 
stem yields an approximate ten percent reduction in the 
flow rate. This equi -percentage valve characteristic is 
obtained by N = 2. This (i.e. Eq, 3*91) essentially simu- 
lates the valve for large valve-stem lifts and pressure 
difference changes. 

3.4.2 PIMP DIRAMIGS 

While simulating the pump trip of the combined 
model of the reactor and turbine, earlier, the term POWP 
in Eqn. 3.11 was set to zero. The neutron power response 
was comparsd with the available test result (Pig .3*2 )• 

Two major differences were observed (a) the dip in the 
neutron power appeared quite early i.e. at less than one 
second, where as in the actual plant it occured at 
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nearly four seconds and (b) the magnitude of the dip l-zas 
much smaller. It was thought tliat this was mainly due to 
ignoring the pump dynamics. Hie pump failure was repre- 
sented as loss of power to the pump (i.e. PO¥P = O). This 
actually means that the pump pressure developed by the pump 
reduces to zero immediately. Tliis however will not happen 
in the actual pump. In an actual pump, flow settles at W/o 
of the value in say 5 seconds depending on pump inertia 
instead of much faster coastdown as in the earlier case. 
This has strong influence on the neutron power variation. 
Also, it is well-established fact that the direct-cycle BWR 
power plant can be controlled through changes in recircula- 
tion pump speed, J0SEFSS}N • Hie above two reasons 
prompted the modelling of the pump dynamics. In boiling 
water reactor the recirculation pumps are driven by induc- 
tion motors. Thus, the pump speed change is a fimction of 
the impressed electrical torque and the mechanical torque 
developed by the pump. Applying Newton' s second law of 
motion to the rotor, the pimip dynamics is represented by 
LINTOED ri6j ■ 


dn 

p 

dt 


30 g. 


TS J. 


P 


(T - T ) 
^ el 


(3.92) 


ytiBTB is tli9 coinbined inoinBnt of inertici of motor^ 

P 

impeller of the pump and other rotating parts# Ihe 
mechanical torque depends upon the pump speed 
the flow rate as shown in Fig# 3 *3 Disregarding 
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small motor time constants, the electrical motor torque is 
obtained from the torque- speed characteristics of the motor. 

-- ■ * ■ ■ Hiese functional 

relationships can be represented by 


T 

m 

- f ^ 

«D- 


(3.93) 

T 

el 

= h 



(3,94) 

graph for 


VS n, 

P 

has not been considered since 


immediately goes to zero with the loss 'of pouer. The 
head developed by the pump itself is a function of the 
pump speed Up and ttie flow rate V/p and is given by 

,.\Pp = f^ (Wp J' (3 *9 5) 

Such a characteristic is given in Fig, 3,4. Thus, the 
Eqn. 3.* 11 can be replaced by Eqn. 3.95. Eor simulation 
purposes, the pump characteristics of Pig, 3.3 and 3.4 are 
resolved into a mesh of T^, A Pp values for flow rate with 
speed as a parameter, Tbese curves are approximated by 
choosing a small enough interval size and the interne- _ 
diate values are calculated by a linear interpolation 
technique. 

3.4,3 SWEEP EFFECTS 

In the reactor sub-system, the core average void 
fraction, and the average quality, Xq, were obtained 
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even during transients by applying quasi-static heat 
balance equations to the core. Thus Eqns. 3.21 and 3.23 
were algebraic in nature and did not take into account the 
dynamic dependenca of voids on flow rate and pressure chan- 
ges. ACKASU was lirst to point out the need, for 

representing the dynamic behaviour. Hiis behaviour has an 
important role to play in transients where large amount of 
flow rate and pressure variations occur e.g., EG pump 
failure. As soon as power to the pump is lost, a rapid 
drop in the flovr rate occurs in the entire primary loop 
including the core. This is accompanied with the reduction 
in heat removal rate from the core and results in the void 
build-up. The neutron power due to reduction in thermalisa 
tion (-ve void coeff. of reactivity), tends to drop. The 
model used earlier did take into account the above features 
However, with the drop in the coolant velocity, the voids 
generated earlier were not being swept as quickly as with 
the full velocity case. Hie sweep phenomenon can be \inder- 
stood with the help of the following arguments. Let us 
watch the moifent of the bubbles in a heated channel shown 
in Fig. 3. 5 . It is assumed that the generation in each 
section is lumped and is the same (generation rate is 
constant) i.e, two bubbles per section Fig, ( 3«5 a) shows 
the void distribution in an initial steady state before 
the pump fails. Now that the velocity drops down to 
half its original value . The sequence of the events at 
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every time step is depicted in Fig. (3.5 b) to Fig. ( 3 . 5 g). 
In constituting this the bubble genemtion rate is assumed 
to hs tile same as earlier i.e, two bubbles per section. 

This is a simplifying assumption and may not be valid in 
practice. In the figures 'O' shows the bubble generated 
at a particular instant of time while shows the bubble 
generated at earlier time but which has moved up. In each 
of the subsequent figures, the bubble marked 1 represents 
the bubble generated at the zeroth step i.e. Fig. (3.5 b). 
It is seen that this bubble marked 1 moves out of the 
channel in six time steps. Further, tiie bubble pattern 
in this channel from the sixth time step onwards is same as 
in fifth time step and the new steady state is reached. 
Since the channel was divided into 3 sections it takes six 
time steps for the bubble marked 1 to escape from the 
channel i.e, it takes the new transit time (calculated on 
the basis of V^/2 and core length) for the bubbles to 
stabilize at the new steady state. In observing the 
Fig. ( 3,5 )j it is clear that the maximum bubble accumula- 
ted is in the 3rd, ^th and 5th time step. This is known 
as tlie sweep effect. 


We will represent the dynamics of the void 
sweep phenomenon by a first order and a single time cons- 


tant dynamics. Thus 
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'^'s 2 

at - Tj 33 ^ ^0 + K (Wjj (t) - Wjj (t - I^g)). CX^) 

(3.96) 

Tug = average transport time across the boiling 

length, secs. 

_ boiling length ^ 

core inlet velocity (3»97) 


The coefficient Kc.' is the partial derivative of void 
fraction o'g, with respect to flow rate, ft^/soc., assumed 
constant and obtainable from a steady-state void map around 
the given operational point. The occurence of Tgg in the 
term (t - ^'p,S^ Eqn. 3 * 9 ^ indicates that the sweep 
effect is felt for a period corresponding to the effective 
passage time for the flow across the boiling length in the 
core. 


The need of sweep effects was strongly felt from 
the study of the results presented in Ref. Hlj, where, 
the void reactivity coefficient was increased five times to 
a pseudo value of 1 to obtain better agreement for 

the pump trip tests. This showed that not all the voids 
have been accounted for in the steady state void repre- 
sentation, whose feedback effect on neutron power is 
most dominant in the present test. Both the 'pseudo* 
void coefficient approach and the explicit modelling of 
sweep have been carried out presently, the latter being 
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seen to give quite meaningful comparisons with the expeii- 
mental resiilts. 


3.4.4- COUPLING OF THE TWO SUB-SYSIMS 


The modified model for the sub-systems i.e., the 
reactor and the turbine are coupled to represent the com- 
plete plant . This is made possible in the model by the 
following features. In the physical process of Fig. 3.1 
the steam drum supplies steam of 99 . 9 ^ quality in the main 
pipe line, a part of which drives the HP turbine and the 
rest is supplied to the shell and tube type reheater for 
superheating the steam coming out of the HP turbine before 
feeding it to the IP and LP turbine. This is represented 
by replacing Eq. (3.36) by 

Wg = C* P (3.99) 


This steam traverses the turbine loop and comes out of 
Heater 2 with certain enthalpy h^^. The Eqn. 3.27 assumed 
a constant feedv/ater temperature which, however, keeps 
varying due to the turbine dynamics as seen in Eqn, 3 . 68 . 
Thus the old Eqn. 3.2? is replaced by the following 
representation 




h ) * 


U 


FW 


^f 


C 

P 


( 3 . 100 ) 


This feedwater, returned from the Heater 2 is fed into 
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st83-in drum vjii©r6 it iiiixes with. 'cIig rscirculo-tion ijatei* wliicli 
is at satrirated conditions. Since these two are at diffe- 
rent temperatures, mixing process in the druia returns the 
mixture at a mixing temperature to the core through down- 
comer. This mixing is not modelled. However, the delay of 
6 seconds associated in the mixing process is modelled as 
a pure transport delay in the simulation. Thus the total 
delay is secs, instead of 9 secs, used earlier. The 
new plant model used for simulation, consists of 21 diffe- 
rential equations and 75 algebraic equations. The differen- 


tial equ£ 

itions a 

re 3.1, 

3.2, 

3.7 to 

3.10, 3 

.28, 

3.3^ 

, 3.^5, 

3.^f6, 3Y 

)2, 3*56 

, 3.57, 

3.60, 

3.61, 

3.6^, 3 

.66, 

3.68 

y3.70, 

3.92 and 

3.96. 
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(3.101) 

and 

u" = 

l.i'^^ex’ 

^el’ 


A 

^2 



(3.102) 


It can be seen that the structures of the equations given 
above are of the type 

X = F jlc, U, Z, X (t - Y )j (3.103) 

where X and U are state and control vectors* 


C3.i<A-) 
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Closer look at the equations would reveal that the terms 
in 2 expressed explicitly in terms of X 8,nd TJj i.e, 

2 = s ( 2, 2 ) 

Therefore substituting eqn. 3. 105 In 3.103 will jdold nonlinear 

differential difference equations in the form of 

X = f ( X , X (t -T ), n) (3.106) 


Simul cl b io n oi the system equations will he discussed in 
next section. 


3,5 SIM'ULATION OF HIE COMPLETE PLANT 

The plant model was simulated by Kunga-Kutta 4th 
order technique, according to GILL The IHi 7044- 

computer took 2.7 secs, for simulating one sec. of the non- 
linear plant model. In order to reduce the computer time 
an approximation was used in the model. This will be dis- 
cussed in Sec. 3«5.2. 

3 . 5.1 VALIDATION OF THE MODEL 

The plant model was initialized to yield a 
corresponding plant operations data before any simulation 
test was conducted. For simulation purposes one of the 
variable operating pressure was selected around 1000 psi, 
since the empirical fits for enthalpy, saturation tempe- 
rature etc, were made in the vicinity of 1000 psi to 
represent Mollier chart and steam table parameters used 
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for calculation. For initialisation the input parameters 
to the reactor e.g. valve lift, motor electrical torque was 
selectoci ]. roni the data while all other parameters were 
calculated from the steady state conditions by equating 
the differential equations to zero. One such operating 
point obtained is given in Table 1. Also the corresponding 
operating data from the plant is shown. The close corres- 
pondence between the two can be observed. The plant model 
is validated only for the Pump Trip, and the Control Rod 
Oscillation tests because of the limited availability of the 
experimental results. Then the following studies are also 
conducted. 


3,5.2 O'fflER STUDIES 

TEEAL OF PROMPT JIMP APPROXIMATION 

VJlien Eq,. 3«1 for the Neutron kinetics is used, 
the integration step size is limited to 0,005 as any larger 
step size would introduce integration error. In order to 
overcome tMs prompt-jump approximation was tried to reduce 
computer time taken for the calculations. The above appro- 
ximation is valid when the relative rate of change of 
reactor power N in Iq, 3»1j ’tbo mean prompt generation 
time 1 , is sufficiently small so that 

1 f m 
^ ^ dt 




1 - 


p 


holds, such that the 
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teim ^ In the kinetic Eqn. 3.1 can be neglected in 

comparison to ( Ak - 1 ) N (t). , This approximation is 
called prompt .approximation because it predicts the 
sudden change in the reactor power following a sudden 
change in reactivity when the concentration of precursors 
remain constant and are represenoed by one group pressure 
model. So Eqn. 3.1 was replaced by linear algebraic 
equation . , 


N = 




(3.107) 


This is described in detail by AICA.SU[62i. This PI appro- 


ximation wan first tried to reactor subsystem using 0,5 


sec, time step, Tiius 100 times faster computer calcula- 
tions wore possible (1 computer sec. = 38 problem sec.), 
introducing less than a maximum of 3% absolute error. 
Figure 3.6 shows the trends for various time steps for 
largo tro.nsiont like pump trip using the PJ approximation. 
But this PJ approximation was not used in the coupled 
system as the integration step size used for the turbine 


part wan the limiting factor as compared to the reactor 
part . 
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TABLE No. 1 

EXPERIMENTAL AND COMPUTED PLANT DATA 


PARAMETER | COMPUTED ' ETJeKIMENTAL 

{ MEASURpM-RNT.q 

1. Neutron Power, Watts • 0,1^7 166 E 09 

2. Delayed Neutrons, " 0.200235 E 11 

3. Average Fuel temp,, °F 0,788573 E 03 

4. COKE Inlet Velocity, ft/sec 0,56833 E 01 

5. COIE imct Tomp. °F 0, 541747 E 03 

6. Reaefeo’c;' . Vesfsol 

ProGsuro psia 0,104384 E 04 

7. Food water weight flow,lbm/^0. 159742 E 03 

8. Steam weight flow, Ibra/sgc. 0,15973 S 03 

9. Recirculation flow rate, 

ft^/sec, 0,772930 E 02 

10, Enthalpy at inlot to HP 

turbine W Sec 0.1254207 E 02 

11, Flow at outlot of HP 

turbine lbs/ sec 0.1423000 E 03 

12, Density of steam in the o 

nozzle Ibs/m-^ O.2325000 E 01 

13, Density of steam at rehea- ' 0.3728000 E 00 

ter exit lbs/m8 

14, Tiithalpy of steam at rehea- 

ter exit W sec, 0.1350000 E 06 

15, Reheating water at reheator 

exit Ib/soc. 0.1149000 E 02 

16, Heat transferred in the 

rdi eater W O.IO19000 E 08 

17 • Flow oi‘ steam at LP turbine 

exit Ib/sec, 0,9200000 E 02 

18, Enthalpy of water at re- 
heator 2 entrance W Sec 0,2357000 E 06 


0.149 E 09 

0.539 E 03 

0.1046 E 04 

0.7816. E 02 


Continued 
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PARAMETER I 

\ COMPUTER 

{EXPERIMENTAL 
MEA RTTRra/IV:: TV .Q 

19. 

Enthalpy of water at 




Reiieater 2 exit W sec 

0.3^5000 E 

06 

20 . 

Cumulative bleed to 
Reheater 2 Ib/sec 

0.3687700 E 

02 

21 . 

Core average exit quality 

0.0593120 E 

00 o,ot 5 

22 . 

Core average exit void 
fraction 

0 . 5 l 2 t 500 

oM 


b, Otiier teats that have been simulated on the model are 
i) 20 '/ and 80?2 step change in the main throttle 
valve, 

ii) valve rod oscillation, 
ill) turbine t'r.f.p . 

3,5.3 CONCLUSION 

Bio tiiodel for the BRP boiling water reactor system 
was presented x-zith a critical approach introducing the 
modifications for better representation of the plant. 

These modifications included throttle valve better repre- 
sentation, the incorporation of pump dynamics and the sweep 
effects. This modified plant was used for simulation 
puriioso. Bio procedure for simulating the pump failure 
tost was (iiscussocl. PJ approximation to minimize the 
computer time was next discussed. 
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FIG.3.2 COMPARISON OF SIMULATION AND PLANT 
TEST RESULTS: FOR RC PUMP FAIL j-r 
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RESULTS, CONCLUSIONS AND RECOI'MENDATIONS 

In this chapter, we present and discuss the 
results of simulation studies for some selected tests 
conducted on BRP boiling water reactor power plant. For 
some of the tests, plant results were available and simu~ 
lation results are compared with those to validate the 
model. Other results for which comparison could not be 
made for lack of test results are also presented as they 
help in understanding the plant behaviour. The chapter 
is concluded with recommendations for further work. 

if.1 STUDIES FOR WHICH EXPERBIENTAL TEST RESULTS V/ERE 
AVAILABLE 

if. 1.1 PUMP TRIP 

It is conducted on the plant by switching off the 
electric power to the pump. This results in the pump speed 
to fall rapidly to the value where it is driven as a turbine. 
This is made possible due to the sufficiently available 
NPSH (net positive suction head) thus the reactor will 
operate on the natural circulation mode. This enables the 
reactor to deliver 50 ^ of the nominal power at the end of 
the transient. 

In the simulation the test was implemented on 
4- versions of model presented in the Chapter 3* They are 
i) without inclusion of the pump dynamics the term 
^Pp (Eqn. 3«11) was set to zero, 


POWP in 
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ii) the pump dynamics was incorporated and the term 
^ei in Eqn. 3.92 was equated to zero, using 
normal void coefficient of reactivity (11.5 jd/fo), 

iii) using pump dynamics and pseudo void-coefficient 
of reactivity (5 times the normal value, refer 
to IBM report h1‘), the term in Eqn. (3.92) 
equated to zero, 

iv) the pump dynamics and sweep model \m.s incorporated 

equating to zero in usual manner, 

for 

Figure 3.2 sliows the response/(A) the neutron 
power, (B) the average heat flux, (G) the recirculation 
flow hy all above mentioned four methods. However, the 
performance for (D) the fuel average temperature and 
(E) the pump speed presented for oiiLy the (iv) case, because 
other methods revealed similar trends. First the plant 
results are discussed followed by the simulation results. 

In Fig, (3,2b) for the neutron power the initial dip is 
observed at nearly 4 seconds then it tries to stabilise at 
75 MW^ after 50 seconds (Fig. 3. 2b), The recirculating 
coolant flow in the core decreases (Fig. 3.2k) result- 
ing in lesser heat removal from the fuel assembly. This, 
essentially, reduces the heat transfer rate Q, although 
at a slower rate. 
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Thus, the average core quality x- increases 

0 

(Eqns. 3.18 and 3.19) increasing the amount of voids in 
the core ( Eqn. 3*21). The increase in voids reduces 

the themalisation action of the moderator in Eqn. 3,1 
resulting in a reduction of Ak the net reactivity of the 
core through (Eq, 3.3)*' This sequence of events conti- 
nues till 6 seconds or so and -the dip in the neutron power 
is noticed. T^ throu^Ja Eqn. 3.7 al^ decreases (Fig, 3,2 p) 
resulting in the reduced neutron ahsorption in the fuel 
due to Doppler effect. The net effect of Akj^p and 
Ak . is that they trj^ to annul the effect of each other, 
till one exactly balances the other and a new steady state 
for the neutron power is reached (Fig, 3*2 ). 

The simulation results also reveal the same nature 
of the test resifLts given, above. The Figtt . 3«2 for the 
neutron power clearly brings out the improvements on 
implementing the four simulation efforts. 

The dip (Fig. 3.2 c) occurs much earlier (1 sec.) 
as compared to 4- seconds in the test result (Fig. 3.2b). 

The incoiporation of the pump dynamics, (Fig. 3 . 2 d) 

(ii) improves the situation with regard to time behaviour, 
(Fig. 3.2 d) while inclusion of the first-order sweep 
effect (iii) gives better magnitude of the dip (Fig, 3 *2 e) 
The fourth method, using pseudo void coefficient, interes- 
tingly, shows an excellent correspondence of the model 
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witli til© plant for 1118 neutron power. This reveals strong 
influence of the voids on the reactor behaviour indicating 
the inadequacy of the first - order sweep. Case (ii) and 
(iii) individually, overestimateswhile (iv) under-estimates 
the time response for' the heat flux. In the actual plant 
also the heat transfer rate is measured using the neutron 
power as an input to the analogue integration scheme 
employing operation amplifiers etc. with a time constant 
of 6.6 seconds as has been done in the model here. Thus 
the above trend for the four different cases can be expl- 
ained due to this. Since tlie reactor pressure was held 
constant with the help of I.P.R. this transient is not felt 
on the turbine side of the plant at all. For other studies 
reported here, the nomal void coefficient of 11.5 Is 
used. Although, the use of pseudo void coefficient exhibi- 
ted excellent correspondence of model with the plant, yet 
there was no strong justification for it to be used for other 
tests. Moreover, during other studies, recirculation 
coolant flow undergoes snaller change and hence the effect 
of void changes due to flow changes is negligible. Thus 
the inclusion of the normal void coefficient does not 
introduce significant errors it was used in other tests. 



h.1,2 CONTROL ROD OSCILLATION 


The control rod oscillation is a steady state test, 
normally conducted on plant to measure the reactivity-to- 
neutron flux transfer fimction by oscillating a standard 
control rod in acertain frequency- range. The steam pressure 
controller, the water level controller and the other control 
rods are kept on manual to avoid variation of pressure and 
terapera.ture due to changes in these control inputs. In 
the simulation this test was conducted by varying the ANkgjj- 
term 2 In Eqn, 3.3 sinusoidally, at a frequency of one 
cycle per sec, (Fig. 1a), llhe neutron power is compared with 
the available BRP plant result, under same conditions of 
magnitude and frequency of the control rod variation. The 
comparison of the simulated and the plant test result for 
the neutron power is shown in Fig. 4-, 1c. 

It is observed that; 

1 . the neutron power of the plant and the . simulation results 
in Figures 4-, 1c and 4.1b varies with same frequency at 
which the control rod oscillates (1 Cps), 

2. the nature of the variation of neutron power (coii 5 )a- 
ring the graphs a and b plotted on same scale) is 
similar for both the^ plant arid the simulation, 

3. There is however a small phase difference of about 

8 - lo'^ between the plant and the simulated results. 
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A niaximimi amplitude of 0.5 MW increase in Q .was observed 
which is negligibly small as compared to 7 MW change in 
the neutron power and also T™ is changed by a maximum 
amplitude of 0.8 (not shown in the figure). This 
perturbation on the plant is hoxirever, damped out before 
it reaches the turbine. 

if. 2 OIHEE STUDIES 

if. 2.1 20^ AND 80^ STEP INSERTION OF VALUE STET-I IN 

THE THROTTLE VALVE 

The throttle waive interconnects the reactor sub- 
assembly with the turbine sub-assembly. In practice, the 
step changes on the valve stem are introduced to study its 
influence on both the reactor and the turbine. 

In the simulation it was performed by setting 
y to 0,20 y^ and 0.80 y^ respectively in the equation 
(3.99) of the following terai which represents the throttle 
valve 

«1 = Sh * % - - h’ 

The responses for these disturbances are shown in Fig- 
if.2 , The comparison with the test result could not 
be done due to the non-availability of the test data. 

The two disturbances of different magnitude carried out 
exhibit similar trends for the plant variables. It is 
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clear from Fig- if *2 that QO% step change has a stronger 
influence on the plant as is to he expected and it is inte- 
resting to study its influence on the plant. As soon as the 
valve stem is dropped.it reduces the area of steam flow. 

The steam flowing out of the steam drum reduces immediately 
and the dip occurs after 0.5 seconds (Fig. 4.2a ), This 
causes the steam flowing in the steam chest to reduce and 
j-Q and h^ as per equations 3.45 and 3.46 suddenly decrease. 
These result in quick fall of the steam chest pressure (Pq) 
as per equation 3.50 with the dip in occuring at 0,5 sec, 
(Fig, 4 , 211 ). The lowering of the steam drawn from the drum 
is felt hy the reactor also which builds up the steam pre- 
ssure through equation 3 * 28 , Thus the saturation ten 5 >eia.- 
ture of the steam increases (equation 3.40) resultirg in 
increase of (equation 3 . 29 ) while (equation 3.18) 

decreases. This, effectively supresses the voids in the 
core (equation 3.21 and 3-23) and the neutron theimalisa- 
tion improves suddenly causing a quick jump in H.at 0.5 sec, 
(Fig. 4. 2d)-. Simultaneously, on the turbine side the 
steam flow to the HP unit decreases resulting in the de- 
crease of hg (equations 3.56 and 3.57) which through 

equations 3.58 and 3.63 decrease Pg and Tg. The response 
for Tg is shown in Fig'- 4.2c ., The steam flow passing 
to the reheater, V/pg shows a small increase due to the 
increase in pressure P. This increased steam flow rate 
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and the decrease in Tj^ mentioned earlier results in an 

increase in as per equation 3.61 (Fig-, lf.2f ). In 

equation 3 *7 "I doth h^ and h 2 increase, the first quantity, becau- 

‘se of linear dependance on pressure P (equation (3.83 )) 

and the latter due to its dependance on (equation 3.^) 

? 

but the h 2 increase is much more compared to hg increase 

I 

so as to decrease (hg “ h 2 ). Also W 2 decreases (not 
shown in Fig, 4.2 ). These result in Tjjp to decrease 

(Fig. 4.2g). The variable hp decreases gradually so 
that T^p (Eqn. 3.42) pursues a gradual decrease (Fig. 

4,2h ), The reactor pressure P keeps building up suppre- 

I 

qding the voids and hence neutron power increasestill the 

, 

average fuel temperature effect is substantial i.e., the 

: ' 'I 

negative feed back effect due to Doppler effect of reac- I 

. ■ li 

tivity. The neutron power continues to build-up slowly till 
15 seconds when the mixing water from the drum flows back 
to the inlet of the core with an Increased sub-cooled 
enthalpy. This causes the voids to increase and hence it 

■ ■ . I 

forces neutron power to decrease. Thus from this stage 
onwards, the net effect on the neutron power will be due 

' * j 

to these three i.e. 

i) flashing effect due to the pressure build-up 
(i.e. void collapse due to the pressure idse),, 

ii) increased void generation rate due to decreased 
subcooling, 

iii) D.oppler effect. 



The net effect results in decreased neutron power. 

The temperature T^^ decreases hecause of the -valv.e 
closure(this is due to the decrease in also), although 
the flow rate passing through the secondary steam valve is 
not changed in this simulation. The decrease in teii 5 >erature 
is observed as shown in Fig,4,2e.Tii8 combined effect of the 
temperature Tj^ decrease and the main steam flow rate decrease 
is ’to increase the heat transfer rate as per equation 
(3»61 ) and is shown in Fig' 4-*2f. Bie respoase of the 
plant is simulated till 20 seconds. It is seen that the 
throttle valve seives effectively to control the turbine 
output. Figt shov/s the decreasing trend of the 

mechanical power produced. 

'4-. 2. 2 TSEOTTLE VALVE BOD OSCILLATION 

Bils steady state test is performed to know the 
variation of the pressure, the neutron power and the 
mechanical power with the variations in the steam flcn^r 
rate . 

In the simulation, it was conducted by sinusoi- 
dally oscillating the valve set point at a frequency of 
1.57 Cps and a magnitude of hofo of rel'^tive position 
through equation 3-91. The responses of the f ollo\'ring : 

A, the relative valve movement, 
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B. the steam flovi rate Wg , 

C. the vessel pressure, 

D. the neutron power, 

E* the heat flux, 

F, the mechanical power are shown in Fig’^ 

The interactions "between various variables has 
already been discussed in Sec. 4,2.1. 

There is an initial transient especially in the 
neutron power Fig-. 4,3d which however dies down in the 
first 3 seconds and afterwards reaches steady state oscilla- 
tion with about 1/2 MW magnitude. While both neutron power 
and reactor pressure Fig, 4.3 are nearly in phase with the 
valve rod oscillation, both mechanical power (Fig,4,3c & e) 
and steam flow rate (Fig. 4.3b ) have - 180 degree phase 
difference with valve rod oscillation. Also a flattening 
of both steam flow rate and mechanical power on the peak 

side of the cycle is observed. This effect for steam flow 

' 2 ' 

is due to the nonlinear valve admittance, (Jq - y ) , 
versus flow characteristics given by Eq. 3.91 and thus 
admittance leads to saturation of the flow rate at high 
flow rates. Since mechanical power directly depends on 
flow rate this flattening is also observed in the mecha- 
nical power curve (Fig', 4.3e) 



B. the steam flow rate Wg , 

C. the vessel pressure, 

D. the neutron power, 

E. the heat flux, 

F. the mechanical power are shown in Fig-^ 

The interactions between various variables has 
already been discussed in Sec* 4-, 2.1. 

There is an initial transient especially in the 
neutron power Fig-. 4-,3d which however dies down in the 
first 3 seconds and afterwards reaches steady state oscilla- 
tion with about 1/2 MW magnitude. While both neutron power 
and reactor pressure Fig, ^.3 are nearly in phase wiiii the 
valve rod oscillation, both mechanical power (Fig, 3c & e) 
and steam flow rate (Fig. 4,3b ) have - 180 degree phase 
difference mth valve rod oscillation. Also a flattening 
of both steam flow rate and mechanical power on the peak 

side of the cycle is observed. This effect for steam flow 

/ \2 

is due to the nonlinear valve admittance, (yg - y ) , 

versus flow characteristics given by Eq. 3,91 and thus 
admittance leads to saturation of the flow rate at high 
flow rates. Since mechanical power directly depends on 
flow rate this flattening is also observed in the mecha- 
nical power curve (Fig'. 4.3e) 



76 


4,2.3 TURBINE TRIP 

Under abnormal conditions occuring on either 
generator, grid or the reactor side of the plant such as 
load thro'w~off etc, the turbine tripping becomes necessary * 
Normally, it is performed by allowing 5^ of the total steam 
through the main steam valve enough to keep the rotor rota- 
ting, and the rest 95"^ being by-passed to the condenser 
through by-pass valve (ref. Fig. 3.1). 

In the simulation, this is performed in tiie follow- 
ing manner. The main throttle valve is closed by 95^ of 
the initial value i.e. y is set equal to 0,95 x y^ to per- 
mit only 5^ of the steam in equation (3.91). In equation 
( 3. 62) representing the secondary valve, the area ±2 is 
reduced to 0,05 ^2 allow only 5^ of the steam flow- 

rate for the purpose of reheating in the reheater. The 
by-pass valve is modelled in such a manner that it allows 
the rest of the steam to the steam condenser initially. 
This is represented by the following equation (4.2) 

^BPV ~ ^BPV ^PV ^ (4.2) 

Where Cgp^ . is evalioated by equating the initial 

95^ steam flow rate to the r.h.s. of equation (4,2) 
at the initial drum steam pressure. This by-pass valve 
is similar to the secondary steam by-pass valve. This 
assimiption is valid only if down-stream pressure i.e. 
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condenser pressure is constant (an assumption already made) 
and the by-pass valve movement does not vary with time 
after the initial step, Wow, in the simulation? eq. (3«99) 
is replaced by the following equation 


W. 


+ Wts-o + W 


'1 


PR 


3PV 


(^.3) 


As is to be expected steam flow from the drum, drops 
down very fast, Fig'. (4.4a) resulting in a steep increase 
of pressure Fig. (4.4 b). Since the steam flow rate VJg 
has decreased, Hie feed water flovj also steeply 
decreases as per Eq, (3*34), Fig , (4.4 c). The increase 

of pressure results in voids to collapse, thus neutron power 
shoots up initially reaching 500 Ml-J (Fig, 4,4 d) at 0,5 
seconds. This causes the heat transfer rate Q and average 
fuel temperature Tj, to increase as sliown in Fig. _ 4,4 e 

and Fig.(4,4 f ) respectively. The increase in T^, brings 
in negative DLoppler reactivity feedback, tlius bringing the 
neutron power back to lower values. Lastly, mechanical power 
Fig.(4,4 g) also decreases steeply since the flow rates 
through HP and LP turbine 'have dscrcashd. 

Some other tests which were simulated includes 
the following 

1, 5/2^ increment in reactivity after 2 seconds, 

2. 10 percent step increase in enthalpy, 

3* pump speed stepped down to 25^j 
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4. pump trip and restart after 10 seconds, 

5. main steam valve closure by 99.8 percent followed 

by two pump trip after 6 seconds, 

they are not reported here as they showed expected beliaviour 
and in order to present the condensed volume. 

4,3 CONCLUSIONS 

A modified mathematical model for the BHP boiling 
water reactor nuclear power plant was studied. The model 
was developed on the basis of the physical laws 
such as mass, energy and momentum balance. This resulted 
in 21 non-linear differential equations and 75 algebraic 
empirical relationships. The model was simulated on the 
digital computer to validate it against the test results. 

The model showed good correspondence with the plant results. 
Some tests for which the plant results were not available, 
were also studied to observe the plant response. The 
emphasis was on studying large order transients on the model, 
however, some small order transients were also studied. This 
study showed that mathematical models could be developed 
to predict the plant behaviour for both small and large 


order disturbances. 
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4-, if EECOMMENDATIONS FOR EXTENSION OP THE NORE 

The work can be extended in two different ways; 

A. SIMULATION STUDIES 

a. Further modifications can be introduced to replace 
the point kinetics by spatial kinetics to study the 
spatial effects* 

b. The voids can be suitably modelled to Include 2nd 
order sweep effects* 

c* The grid generator sub~syste]ii can be coupled to 

the above plant to study the overall plant response. 

d. Variable time delays especially for pump trip tests 
etc. sho-uld be modelled. 

B. MODERN CONTROL THEORY 

a. The plant model can be linearised and the feed back 
control approach can be applied to design the 
feedback controllers by using optimal modal control 
techniques. 
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